Global regularity and stability of a hydrodynamic system 
modeling vesicle and fluid interactions 



Hao Wu * and XiANG Xu ^ 
February 23, 2012 

(N 
O 

<n: 

^ , Abstract 

I I In this paper, we study a hydrodynamical system modeling the deformation of vesicle 

membranes in incompressible viscous fluids. In the three dimensional case, we prove the 
P\j . existence/uniqueness of local strong solutions for arbitrary initial data as well as global 

strong solutions under the large viscosity assumption. We also establish some regularity 
P I ■ criteria in terms of the velocity for local smooth solutions. Finally, we study the stability 

■ of the system near local minimizers of the elastic bending energy. 

|: 

1 Introduction 

Biological vesicle membranes are interesting subjects widely studied in biology, biophysics 
^ ■ and bioengineering. They are not only essential to the function of cells but exhibit rich 

physical and rheological properties as well [23]. The single component vesicles are possibly 



OO ■ the simplest models for the biological cells and molecules, which are formed by certain amphi- 

■ philic molecules assembled in water to build bi-layers [12] . The equilibrium configurations 

■ of vesicle membranes can be characterized by the Helfrich bending elasticity energy of the 
CN ■ surface [HllTlfTB] such that they are minimizers of the bending energy under possible constraints 

like prescribed surface area and bulk volume that account for the effects of density change and 
osmotic pressure [T2|l32] . Let T be a smooth, compact surface without boundary representing 
the membrane of the vesicle. In the isotropic case, if the evolution of the vesicle membrane 
does not change its topology, the interfacial energy takes the simplified form [7]: 



X 



-^elastic - j^ij^i^ - Hof ds, 



where H is the mean curvature of the membrane surface; k is called the bending rigidity, 
which can depend on the local heterogeneous concentration of the species; Hq is the spon- 
taneous curvature that describes certain physical/chemical difference between the inside and 
the outside of the membrane. 
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Recently, phase-field models have been derived within a general energetic variational 
framework to study vesicle deformations and numerical simulations of the membrane defor- 
mations were carried out (see e.g., [Tm - [T5l[3n[32] and references cited therein). As in [91IT2]. 
we denote by (j) the phase function defined on the physical domain 0, which is used to la- 
bel the inside and the outside of the vesicle F such that (p takes the value 1 inside of the 
vesicle membrane and —1 outside. The sharp transition layer of the phase function gives 
a diffusive interface description of the vesicle membrane F, which is recovered by the level 
set {x : <^(x) = 0}. The phase field approach describes geometric deformations in Eulerian 
coordinates and it provides a convenient way to capture topological transitions such as vessel 
fission or fusion via changes in the level set topology. This simplifies numerical approxi- 
mations because it suffices to consider a fixed computational grid rather than tracking the 
position of the interface 

For the sake of simplicity, we assume that A; is a positive constant and Hq = 0. The phase- 
field approximation of the Helfrich bending elasticity energy is given by a modified Willmore 
energy |12^l32j (see, e.g., [TO] the approximation energy for the elastic bending energy with 
non-zero spontaneous curvature) 



where e is a small (compared to the vesicle size) positive parameter that characterizes the 
transition layer of the phase function. The convergence of the phase-field model to the original 
sharp interface model as the transition width of the diffuse interface e — > was carried out 
in [9l[32]. Two constraints are widely used in the biophysical studies of vesicles [26] such 
that the total surface area and the volume of the vesicle are conserved (in time). The former 
is a consequence of the incompressibility of the membrane, while the latter is based on the 
consideration that, for a fluctuating vesicle with the inside pressure and outside pressure 
balanced by the osmotic pressure, the change in volume is normally a much slower process 
in comparison with the shape change [11] . The constraint functionals for the vesicle volume 
and surface area are given by (cf. [T2] ) 



Two penalty terms are introduced in order to enforce these constraints, and the approximate 
elastic bending energy is formulated in the following form [9lll0 1 [T3 l ll4]: 



where Mi and M2 are two penalty constants, a = A[(j)Q) and /? = B{(j)Q) are determined by 
the initial value of the phase function 4>q. Alternatively, Lagrange multipliers could be used 
to conserve the vesicle volume and total vesicle surface area [TTlfH] . 

In this paper, we consider a hydrodynamic system for the interaction of a vesicle with the 
fluid field, which describes the evolution of vesicles immersed in an incompressible, Newtonian 
fluid [8]. More precisely, we study the following phase-field Navier-Stokes equations for the 
velocity field u of the fluid and the phase function (p: 







E{^) = EM + \mi{A{4>) - af + \M2{B{(t>) - (3f 



(1.3) 



ut + u-Vu + VP = nAu + 



54) 
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V ■u = 



(1.5) 



(pt + u- V(/> = -7 



5(1) 



(1.6) 



System (jl.4p ~ ()1.6p can be obtained via an energetic variation approach [M] (see [TT] 
for the derivation of a corresponding evolution system that adopts the Lagrange multipher 
approach for the volume and surface area constants). The resulting membrane configuration 
and the flow field reflect the competition and the coupling of the kinetic energy and membrane 
elastic energies. Equations ()1.4p and (jl.Sp are the Navier-Stokes equations of the viscous 
incompressible fluid with unit density and a force, which is derived from the variation of 
the elastic bending energy and it involves a nonlinear combination of higher-order spatial 
derivatives of the phase function. P denotes the pressure and ^ is the fluid viscosity, which 
is assumed to be a positive constant throughout both fluid phases and the interface. ()1.6p is 
a relaxed transport equation of (j) under the velocity field u. Its right-hand side contains a 
regularization term, where 7 is the mobility coefficient that is assumed to be a small positive 
constant. 

Well-posedness of the system (jl.4p ~ (jl.6p subject to no-slip boundary condition for the 
velocity field and Dirichlet boundary conditions for the phase function has been studied 
in [8l[22]. In [8], the authors obtained the existence of weak solutions by using the Galerkin 
method. They also obtain the uniqueness of solutions in a more regular class than the one 
used for existence. Quite recently, existence and uniqueness of local strong solutions have 
been proved in [22] via a fixed point argument. Existence of almost global strong solutions 
is obtained under the assumption that the initial data and the quantity + a)^ are small 
enough. However, they were not able to prove global existence result because uniform-in- 
time a priori estimates were not available in their argument. On the other hand, since some 
compatibility conditions (at the boundary) are required in the fixed point strategy described 
in [22] to obtain enough regular solutions, the authors have to confine them to function spaces 
with proper limited regularity. 

We note that, although Dirichlet type boundary conditions are more natural and physical 
conditions, the periodic boundary conditions can also be reasonably justified physically when 
the vesicle interface T is sufficient small compared with the overall physical domain Q, (cf. |14)). 
In our present paper, we study the system (jl.4p - (jl.6p subject to the periodic boundary 
conditions (i.e., in torus T^): 



where Q is a unit square in M^. 

The main purpose of this paper is to study the existence, regularity and stability of global 
strong solutions to problem (jl.4p - (|1.8p . In the subsequent proof, we shall see that problem 
(jl.4p - (jl.8p has an energy dissipation mechanism (cf. (j2.16p below) that plays a crucial role 
in controlling the contribution to the momentum equation of the extra stress tensor due to 
the membrane deformation and the contribution of the convection term to the phase-field 






and to the initial conditions 



u\t=o = uo{x), with V • Mo = 0, / uodx = and (j)\t=o = (poix), for x £ Q, 

Jq 
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evolution. The advantage to work in the periodic setting is that one can get rid of certain 
boundary terms when performing integration by parts. Due to the weak couphng in the 
phase-field equation (|1.6p that is a gradient flow of the elastic bending energy under the 
fluid transport, we can derive uniform-in-time estimate for i/^-norm of (j) (cf. Proposition 
13. 1|) that enables us to derive some specific higher-order energy inequalities (cf. Lemma 13.31 
and Lemma l3.4p in the sprit of [21] for a simplified nematic crystal system. Based on these 
higher-order inequalities, we can show existence and uniqueness of local strong solutions to 
problem (|1.4|) ~ p.8|) (cf. Theorem 13. ip . existence of global strong solutions under properly 
large viscosity (cf. Theorem 13. 2p and also the eventual regularity of the global weak 
solution (cf. Corollarv 15. ip . After a careful exploration of the nonlinear coupling between 
velocity field and membrane deformation, we establish some regularity criteria for solutions 
to problem (jl.4p ~ (jl.8p in 3D that only involve the velocity field (cf . Theorems 14.11 14. 2p , 
which coincide with the results for conventional Navier-Stokes equations. This indicates that 
the velocity field indeed plays a dominant role in studying regularity for solutions {u,(j)). 
Finally, we prove the well-posedness and stability of global strong solutions (cf. Theorem 
15. ip when the initial datum is close to a certain local minimizer of the elastic energy by 
using a suitable Lojasiewicz-Simon type inequality (cf. Lemma l5.2p . The results obtained 
in this paper hold for any given (but fixed) penalty constants. Since now we are working 
with the penalty formulation to incorporate the volume and surface area conservation of the 
vesicle membrane, the constraints are satisfied only approximately. It would be interesting 
to investigate the corresponding results for the evolution system in the Lagrange multiplier 
formulation (cf. [11]) where the constraints are satisfied exactly. 

The rest of the paper is organized as follows. In Section 2, we present the functional 
settings and some preliminary results. In Section 3, we prove the existence of local strong 
solutions and global ones under the large viscosity assumption. In Section 4, we establish 
some logarithmic-type regularity criteria for the smooth solutions only in terms of the velocity 
field. In Section 5, we study the well-posedness and stability of global strong solutions near 
local minimizers of the elastic energy. In the final Section 6, we sketch the proof of the 
Lojasiewicz-Simon type inequality that plays a key role in the analysis of Section 5. 

2 Preliminaries 

We recall the well established functional settings for periodic problems (cf. [29]): 

{u £ H"'(W\R) I u{x + ei) = u{x)}, 
{v £ L2(Q), \7.v = 0}, where ^(Q) = H°(Q), 
{veuliQ), \/-v = 0}, 
the dual space of V. 

For any Banach space B, we denote by B the vector space {By, r € N, endowed with the 
product norms. For any norm space X, its subspace that consists of functions in X with 
zero-mean will be denoted by X such that X = (z X : JqW dx = 0^ . We denote the inner 
product on Lp{Q) (or Lip{Q)) as well as H by {■,■) and the associated norm by |[ • ||. The space 
HTiQ) will be short-handed by HIJ^ . We denote by C and Cj, i = 0, 1, • • • genetic constants 



H = 
V = 

v = 



4 



which may depend only on fi,j,Q,a,f3 and the initial data (no,(/>o)- Special dependence 
will be pointed out explicitly in the text if necessary. Throughout the paper, the Einstein 
summation convention will be used. Following [29], one can define mapping S 



Su 



-Au, Vu G := {u G H,Au G H} = nH. 



(2.9) 



The Stokes operator S can be viewed as an unbounded positive linear self-adjoint operator on 
H. If D{S) is endowed with the norm induced by L^, then S becomes an isomorphism from 
D{S) onto H. More detailed properties of operator S can be found in [29j. We also recall the 
interior elliptic estimate, which states that for bounded domains Ui CC U2 there is a constant 
C > depending only on Ui and U2 such that ||</'||ff2(f/j) < C{\\A(f>\\l2(^^^) + |i</'||L2([/2)). In 
our current case under periodic boundary conditions, we can choose Q' to be the union of Q 
and its neighborhood copies. Then we have 



miQ) < C(||A(^|U2(Q,) + Ml^q')) = 9C{\\Acj,\\LHQ) + Ml^q))- 



(2.10) 



It follows from the periodic boundary condition that Jq V(j)dx = and Jg A(j)dx = 0, then 
we infer from the Poincare-Wirtinger inequality that 



IIV0II + 



H^, l|A(/'|| + 





>dx 


Jq 





H2, \\VAcj)\\ + 



(pdx 



HZ. (2.11) 



A direction calculation yields that the variation of the approximate elastic energy is given by 
5E{(t)) 



kg{<t>) + Mi(A(0) - a) + M2{B{^) - (3)f{^), 



(2.12) 



where 



5(</.) = -A/(,^) + l(3</.2-l)/(</>). 



Since we are now dealing with the periodic boundary conditions, the average of velocity u is 
conserved. 



Lemma 2.1. Let {u,(f)) be a solution to the problem ()1.4p - ()1.8p on [0,T]. It holds 



u{t)dx = / u^dx, VtG[0, T]. 
Q Jq 



(2.13) 



Proof. It follows from ([232]) that 

= kg{(P)V(j) + Mi{A{(P) - a)V(j) + M2(5(0) - /3)/(</>)V0 := h + h + h- (2.14) 



Since A{<j)) and B((p) are functions only depending on time, using integration by parts and 
the periodic boundary conditions, we deduce that 



I2 = Mi{A{^) - a) / V^dx = 



I3 = M2{B{<P) - f3) / f{cP)V(l)dx 
Jq 



5 



M2(i?(0) - /5) V (||V</.|2 + ^(0' - 1)') dx 

-M2iB{(j)) - f3) [ V • (V0 (g) V0)dx 
JQ 



= 0, 

where we have used the fact A(pV(j) = V 
1 



k 



h 



£ / A^4>V(t>dx - - 
Q £ 



+ 



1 



V(f)) - iV(|V(/)|2). Finally, 

A(03 - (t))V(l)dx 
\ J AcpV{cP^-cl))dx + ^ J V[{cp' - cPf]dx 



e 



-eA0+ - 1)0 



J / V\Acp\^dx-- 
2 . n £ 



0. 



Thus, we conclude that 



-— Vffl dx 



0. 



(2.15) 



After integrating (|1.4|) over Q, we infer from (|1.5|) . the periodic boundary condition p.7|) and 
(ITT^ that (imjl holds. □ 

Remark 2.1. By Lemma \2.1\ if one assumes that the average of the initial velocity vanishes, 
i.e., 1^ fgUodx = 0, then we can apply the Poincare-Wirtinger inequality to the solution u 
such that the H^-norm of u can be controlled by ||Vn||. When a flow with non-vanishing 
average velocity u is considered, as for the single Navier-Stokes equation (cf. JM3i), we can 
introduce the variable u = u — Jq udx and transform the problem p.4p - (|1.8p into a new 
one in terms of u and (f). Since J^udx is a known constant determined by (j2.13p . it is 
not difficult to verify that our results on existence and uniqueness of weak/strong solutions 
for the initial velocity with zero mean can be extended to this case with minor modifications. 
However, results on long-time dynamics in Section 5 are no longer valid, because the velocity 
will not decay to zero (we also refer to I33j for a similar situation for the liquid crystal 
system). 

For the sake of simplicity, in the remaining part of this paper, we will assume that the 
average flow vanishes. An important property of the coupling system ()1.4p - ()1.8p is that it 
has a basic energy law, which indicates the dissipative nature of the system. It states that the 
total sum of the kinetic and elastic energy is dissipated due to viscosity and other possible 
regularization/relaxations rates. A formal derivation can be carried out by multiplying (|1.4p 
by n, (|1.6p by ^t^, respectively, and integrating over Q. As a consequence, we have 
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Lemma 2.2 (Basic energy law). Let {u,4>) be a smooth solution to the problem (|1.4p - ()1.8p . 
The following dissipative energy inequality holds: 



d /I 



dt \2 



\u{t)f + E{(P{t))^ +^||Vnf + 7 



5E 



0, Vt>0. 



(2.16) 



Based on Lemma 12.21 we can apply the Galerkin method similar to that in [8] to prove 
the following result on existence and uniqueness of weak solutions to the problem ()1.4p - ()1.8p . 

Theorem 2.1 (Existence of weak solutions). Let n = 3. For any initial datum (mqi'/'o) ^ 
H X T > 0, there exists at least one global weak solution {u, cj)) to the problem ()1.4p - (jl.8p 
that satisfies 



u G L°°(0,r;H) nL2(0,T;V); 

€ L^{0,T;H^)nL\0,T;H^)nH\0,T;Ll) 

In addition, the weak solution is unique provided that u G L^{0,T;'L'i). 



(2.17) 
(2.18) 



Besides, we can obtain the following uniform-in-time estimates on weak solutions from 
the basic energy law: 

Proposition 2.1. Suppose n = 3. For any initial data uq G il, (po £ Hp, the corresponding 
weak solutions of the problem ()1.4p - ()1.8p have the following uniform estimates 



\um + \m)\\H^<c, vt>o, 



(2.19) 



fi\\Vu{t)f + -f 



5E , , 



dt < C, 



(2.20) 



where C > is a constant depending on ||tio||, ||0o||_ff2 o-nd coefficients of the system except 
the viscosity fi. 

Proof. We can derive a weaker version of the basic energy law rigorously via the Galerkin 
procedure such that the weak solution (n, (p) to problem (|1.4p - (|1.8p satisfies 



1„ 
— n 



+ Em))+ / ^||Vn(s)f + 7 



5E_ 



ds < ^\\uof + E{ 



yt > 0. 



Recalling the definition of E, we know |||tio|P + E{0) can be estimated by a constant de- 
pending on ||uo||, |j</'o||_f/2 and coefficients of the system, but not on /x. Thus (|2.20p holds and 
||u(t)||, E{t) are bounded. On the other hand, we infer from the boundedness of E(t) that 



l^((^)l<c, \B{<p)\<c, \\fm<c, yt>o. 

Hence, | Jq (j)dx\, \\V4>\\ are bounded. Then by the definition of f{(j)) and Sobolev embedding 
theorems, we can deduce that ||(/'||_h-2 is bounded. The proof is complete. □ 
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3 Existence of strong solutions 



In this section, we study the existence of strong solutions. For this purpose, it suffices to derive 
proper higher-order uniform estimates for the Garlerkin approximation of weak solutions and 
then pass to limit. We observe that the entire calculation is identical to that as we work with 
classical (smooth) solutions to problem ()1.4p - ()1.8p . Thus, for the sake of simplicity, all the 
calculations below will be carried out formally for smooth solutions. 

By the Sobolev embedding theorem in 3D, we can derive the following estimates that will 
be frequently used later. 



Lemma 3.1. Suppose n = 3. We have 

6E 



IIVA0II < C 



2 „ 9 „ 1 

+ C, ||A2,^i|< — 



5E 



+ C, V(AGi/;, 



where C is a constant depending on ||(/>||//2 and coefficients of the system. Besides, 



IVA^^II < ± 



6(t> 



+ C, 



where C is a constant depending on \\(j)\\H3 and coefficients of the system. 
Proof. Recalling ()2.12p . we can rewrite — ^.s 



5E{(I)) 



(3.1) 



where 



+ ^(30 



--</<|V<AP--(3<A^-1)A0 + A(3<^2 
e e E'^ 



2 l)/(0) + Mi(^(</)) -a)+ M2{B{<P) - I3)fi<p) 



+Mi(^(</>) - a) + M2(i?(0) - /3)/(</>). 



(3.2) 



By the Holder inequality and Sobolev embedding theorems, we infer that 



A" 



< — 



C||</.||i.o||A,/.|| 



+ c\\^L^\\v^\ 

+C||A0|| +C(||0||ioo + 1) + Afi||</.||ii +Mia 

+cM2(/3 + \\v^^ + c\\m, + c)i\\A^ + cuifie + cm 



< — 



1 

ke 



6E 



(3.3) 



where C is a constant depending on ||0||j:f2 and coefficients of the system. The estimate 
for \\VA(j)\\ easily follows from 1^ and the fact ||VA0f = \J^A(j)A'^(l)dx\ < C\\A'^(j)\\. 

2.1, ■ , ^ ^E{(P) , ^ . 

I ^^11 -^g j^gt apply V to — — — and we can obtain our 



Concerning the estimate for ||VA 
result by estimating 'SJH{(j)) via proper Sobolev embeddings. 



□ 
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Since our system ()1.4p ~ ()1.6p contains the Navier-Stokes equations as a subsystem, in the 
three dimensional case, one cannot expect that the weak solutions will become regular for 
strictly positive time. But it is worth noting that, due to the weak coupling in the phase-field 
equation ()1.6p that u enters in the evolution equation only as a lower order term u ■ Vcp, we 
can first derive certain regularity results for the phase function <j) and show that it turns out 
to be regular for t > 0. 

Lemma 3.2. Let n = 3. For any smooth solution to the problem ()1.4p - (jl.Sp . it holds 

^||VA(/>f + kJe\\VA^^^ < C{\\Vuf + l)|lVA</.f + C(l + \\Vuf), (3.4) 

where C > is a constant depending on ||tio|l, ||(/'o 11/^2 and coefficients of the system. 
Proof. Multiplying ()1.6p by —A^cp, integrating over Q, we have 

^|||VA0f = -(V(n . V<A), VAV) +7(V^, VAV). (3.5) 

Using the lower-order uniform estimates in Proposition [JTll we estimate the first term on the 
right-hand side of (|3.5p as follows 



(V(n- V</.),VA 
< ^llvA^f + C\\Vu ■ V0f + C\\u ■ V2^"2 



< ^\\vA^^^ + c\\vuf\\v^l^ + c\\u\\l4v^^l, 

o 

< ^!|VA2</>||2 + C||Vnf ||V<^||hi||V</>||h2 

< ^IIVAVIP + C\\Vufi\\VA(l)f + 1). 



For the second term, we infer from (|2.12p that 
7fv^,VA^ 



7 (v(%(0) + Mi{A{ct>) - a) + M2{B{4>) - P)f{4>)) , VA^) 
kjiS/Af{^),VA^cP) - ^(V[(302 _ vA^c/,) 



where 



+M27(i?(0)-/3)(V/(0),VA2<; 
h+h + h, 



h < -fc7e||VAVf + y||VA((/.3-0)||||VA20|| 

< -I^llVA^^f + C\\VA^' + Cll^lli. II A^llia 
+C||</-||ioo||VA0f + C||V0||66 

< -^l|VA20||2 + C(||VA0f + l), 
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h < 



< 



< 



kje 



{\\Gfi^)cpVcP\\\\VA^<P\\ + 11(302 - 1)V/((A)||||VA^ 

vA'^'+c\\A^ls\m\leml^+c\\v^' 



+cw-l\\l^{\\vA^^ + \\cp^- 

^||VA20||2 + C7(||VA0f + l), 



|2 

Il° 



IIVcAf) 



/3 < ^llVA^^f + C(i?(,/.)-/3)2||V/(</.)f 

< :^||VA2</,||2 + C(|lV</)r + 110' - 111' + 1)(||VA0||2 + 1102 _ l||2^ ||V0||^ 

< ^||VAV||2 + C7(||VA0f + 1). 

o 

Collecting the above estimates together, we arrive at our conclusion ()3.4p . 



□ 



Based on the higher-order differential inequality (|3.4p for 0, we get 



Proposition 3.1. Suppose n = 3. For any uq € H, 0o € Hp, the weak solution to problem 
(fLiD^ ffLHjl satisfies 



< C 1 + 



t 



and ||V0(t)||L- < C ( 1 + - ) , Vi > 0, 



(3.6) 



where C is a constant depending on \\uo\\, ||0o 11/^-2 and coefficients of the system. Moreover, 
if we assume in addition that 0o € Hp, then 

||0(t)lb3<C7 and ||V0(t)||Loc < C, \ft > 0, (3.7) 

where C is a constant depending on \\uq\\, ||0o||h3 and coefficients of the system. 

Proof. We infer from Proposition 12.11 and Lemma l3. II that for any r > and t >0, 



j-t+r rt 

sup / ||VA0(r)f dr < supC / 

i>0 Jt t>0 Jt 



rt+r 


SE 




Jt 





dr + Cr 



< C 



r+oo 


6E 






Jo 





dr + Cr < C(l + r) 



rt+r r+oo 
sup/ \\Vu{T)fdT < \\Vu{T)fdT <C. 

t>0 Jt Jo 



\\/u{T)\\^dT < 

Then it follows from (13.41) and the uniform Gronwall lemma 1301 Lemma III.l.ll that 



||VA0(t + r)f <C(1 + ^), Vt>0, 



(3.8) 
(3.9) 

(3.10) 



which yields (|3.6|) . The estimate for ||V0(i)||L°o follows from the continuous embedding 

If we further assume that 0o G Hp , then by the standard Gronwall inequality, we see that 
||VA0(t)|| is also bounded for t € [0,1]. This combined with (|3.10p yields our conclusion. 
The proof is complete. □ 
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Remark 3.1. We remark that the generic constant C throughout the proof of Lemma \ 3.S\ 
does not depend on the viscosity fi, thus the uniform bounds for ||</>||//3 obtained in Proposition 
3.1\ is independent of /x. 



Define 



A{t) = \\Vuf{t) + 7] 



5E 



where r/ > is a proper constant to be determined later, which might depend on H-uqI 
and coefficients of the system. 



(3.11) 

4>o\\m 



Lemma 3.3. Let n = 3. For any smooth solution to the problem 

110(011^3 + ||V0(t)||L- <i^, Vt>0, 

then for 



^f 



V 



^7 



the following higher-order energy inequality holds: 



A— 



<CM^{t)+A{t)), 



(3.12) 
(3.13) 

(3.14) 



where C=k is a constant depending on \\uq\\, \\(j)Q\\ij2, K and coefficients of the system. 
Proof. By equation (|1.4p and the periodic boundary conditions (jl.7p . we can see that 



— r Vti 



-{ut, Aii) = — /x|| Ati|p + (n • Vn, An) 



^5E 



V(/), An 



(3.15) 



Using the uniform estimates ()2.19p and ()3.12p . the right-hand side of ()3.15p can be estimated 
as fohows 

(n-Vn,An) < — IIAnlP + (7||n • VnlP < — ||An|P + C'||n||r oc ilVnll^ 

< — IIAnf + C7(||Vn||||An|| + ||Vnf)||Vn||2 
16 

< ^||An||2 + C(||Vnf + llVnf), 



(3.16) 



r^V,^,An) < ^||Anf + ^ 



5E 



WVct^Wioo < ^IIAnf + C7 



\0<p J » ■ ■ /X 

On the other hand, using integration by parts, we obtain from ()1.6p that 

bE 2 



bE 



/6E d 6E\ 



ld_ 

, /a2. /c/ 5E\ Qk f , 5E 

fe(AV.^)--(9.[A(^^-<^)],^j+^ 



4((3^^-W(^),f)+M,|A(^) 



5E 



dx 



+M2j^B{4>) (f{^), ^) + M2(i?(0) - /3) (dtf{^), ^ 



(3.17) 
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1=1 

The first term Ji can be estimated as follows 



(3.18) 



< 



< 



< 



-kej 

7ke-/ 



6E 
A— 



< - 



Ykej 



< - 



Ike-f 



7fc£7 



5E 
1^ 
5E 
1^ 
5E 

A— 
1^ 



-ke (A(n-V<^),A^ 

V^||A(..V0)f 

+ —{\\Au ■ V(/.f + 2||Vn ■ V'^^f + \\u ■ VA.^f ) 

+ _(||Auf ||V<^||i^ +2||V<3||</)f^,,e + ||K||i^||VA<^f ) 
2 2hFK'^ 

+ — — ||A«f + CK2||A^x||||V^^|| 

7 

+ ^^ff^||Auf + C||Vu|p. (3.19) 

7 



Then for J2, J3, ^4, a dhect computation yields that 



^2 + ^3 + J A 

6k 



6E\ 3k 



= dti\V(l>\'cP), -- - - dt{(t>'Acj^), -- + - ^<Pt, TT 



6E\ k 



6E 



6k / ^ 6E\ 6k /,,4 
--^^A#„-j+-3^(<^^ 



e 

6k 

£ 



6E\ k 



(3</>^ - l)A0t, ^ + ^ (3'/' - 1) 



12A: 



6E\ 



5E\ 



+^ ([2(1 - 3^')A<f>t - 6V<^2 ■ V</>t], ^) 



«^2a + «^26 + <^2c + <^2d- 



(5^ 



Then we have 



J2a < C||V</.||ioc||</.t 

5E 2 



5E 



< C\\u-V(f)\\ 



6E 



5(1) 



+ C 



6E 



< C 

< C 



L3 



6E 



+ C7||Vn|^ 



J2b < C||</>||z,^||<^t||||A</>||^6 ^ 

/ ^fu II iiw^ii SE \ f ^5E \ 5E I SE \ 
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< 



A 



5E 



+ C 



5E 



C\\Vu\ 



J'. 



2c 



< 



< cml^ + m^4>t\\ 



< c 



< c 



A— 



5E 



6E 



ke'f 



A 



6E 



+ C 



+ 

5E 2 



+ C\\V{u-V(l))f 

+ c\\vuf\\vm^+c\\u\\l4 

+ C\\Vuf, 



< cm\u + i)Ut\\ 



5E 



6(t> 



< C 



5E 



+ C\\Vu\ 



Hence, we obtain 



J2 + J3 + J4<—r 



A— 



For the remaining terms, we have 



J^<c\\4>t\\ 



6E 



< C 



+ C 



6E 



6E 



+ C\\Vuf. 



+ C\\S/uf, 



Je < C(||V</.||||V</.i|| + ||(/>3 



6E 



< 



kej 



6(p 



+ C 



5E 



+ C\\Vuf 



Jr = -M2e(i?(0)-/3)(A<Ai,^) + ^(i?(0)-/3)((3</.2-l),/,,,^) 



ke^ 



< 



5E 



A— 



+ cUt\\ 

5E 2 



5E\ Ah 

e 

6E 



M2e{B{4>) - P) A^) + -^{B{<t>) - /3) ((S^^ _ ^^^^^ 6_ 
5E 



5E^ 



+ C 



6(t> 



+ C\\Vu[ 



Collecting the above estimates, we deduce that 



1 d 
2di 



5E 2 ke-f 



2 AkeK"^ 



< 



7 



\\Auf + C ^ +C\\Vu\\ 



Now multiplying (j3.20p by rj 
The proof is complete. 



fi-y 



(3.20) 



and adding the result to ()3.15p . we obtain (j3.14p . 

□ 



13 



Theorem 3.1 (Local strong solution). Let n = 3. For any initial datum (uo,0o) S V x 
Hp{Q), there exists Tq € (0, +00) such that problem ()1.4p - ()1.8p admits a unique strong 
solution {u, (j)) satisfying 

uGL°°(0,ro;V)nL2(0,ro;H2); (3.21) 
^ G L^{0,To;H^)nL\0,n;H^)nH\0,To;H^p). (3.22) 

Proof. It follows from Proposition l3.1l tliat the assumption (|3.12p in Lemma r3.3l is satisfied and 
K is a constant depending on ||tio||, l!'/'o||_ff3 and coefficients of the system. As a consequence, 
(|3.14p holds with t] and C* depending on ||no||, ||0o||_ff3 and coefficients of the system. A 
standard argument in ODE theory yields that there exists a Tq = ro(^(0),C*) G (0, +00) 
such that A{t) is bounded on [0, Tq]. The bound only depends on Tq, A{0) and C*. This fact 
together with the lower-order estimates in Proposition 12.11 and the Galerkin scheme similar 
to that in implies the existence of a local strong solution to problem p.4p - p.8p in the 
time interval [0,to]. Since u G L°°(0,ro;V) C L^{0,To;^), uniqueness of the local strong 
solution follows from Theorem 12.11 The proof is complete. □ 

In general, we cannot expect existence of global strong solutions to problem ()1.4p - (|1.8p 
for arbitrary initial data in V x Hp, due to the difficulty from the convection term in the 3-D 
Navier-Stokes equations. However, if we assume that the viscosity fj, is properly large, then 
problem ()1.4p - ()1.8p will admit a unique global strong solution that is uniformly bounded 
in X H'^ on [0,+cx3). To verify this point, we first derive an alternative higher-order 
differential inequality. 

Lemma 3.4. Let n = 3. For arbitrary /ig > 0, if > iJ-o > 0, and ()3.12p is satisfied, then 
choosing the parameter r] in A{t) to be 

the following inequality holds for the smooth solution {u,<j)) to problem ()1.4p - ()1.8p 



d 



5E 



2 



A{t) + ( n - n2 A{t)] \\Auf + ke-/r]' " < CU(t), (3.24) 

dt \ / d(p 

where C is a constant depending on \\uo\\, ||(^o||_ff2, K, ^0 o.'ri'd coefficients of the system but 
except fi. 

Proof. We only need to refine the estimate (|3.16p in the proof of Lemma 13. 3[ 

(n-Vn,An) < -HAuH^ + -HuHloo ||Vn|p 
8 jji 

< -\\Auf + — (||u||^||Ati||t + )||Vuf 
8 fi 

1 

< ^llAnf + ^llVnf IIAnf + C (/.-^ + fi'') \\\/uf. (3.25) 
8 2 V / 

Since n > Ho, we can choose rj in A{t) to be rj' = j^^j^- Combing (j3.25p with estimates 
for the other terms in the proof of Lemma 13.31 can easily conclude (|3.24p with our choice 
ry'. □ 
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Theorem 3.2 (Global strong solution under large viscosity). Let n = 3. For any initial data 
{uo,do) G V X Hp, if fi is sufficiently large such that ()3.26p is satisfied, then the problem 
()1.4p ~ (|1.8p admits a unique global solution. 

Proof. We infer from ()2.20p and the choice of ij' in Lemma 13.41 that 

rt+l r+oo 
sup / A{T)dT < / A{t)dt < M, 

t>0 Jt Jo 

where M is a constant depending on ||uo||, ||<;/'o||/i'3, fJ-o and coefficients of the system but 
except /J,. If the viscosity fi satisfies the following relation 

> A{0) + Cm + 4M + /i| , (3.26) 

by applying the same idea as in [211 Section 4], we can deduce from p.24p that A{t) is 
uniformly bounded such that 

A{t)<fi^, Vt>0. 

Based on the uniform-in-time estimates and the Galerkin scheme, we are able to prove the 
existence and uniqueness of a global strong solution to problem (|1.4p - ()1.8p . We leave the 
details to interested readers. □ 



4 Regularity criteria 

In this section, we are going to establish some regularity criteria for solutions to problem 
(jl.4p - ()1.8p in the three dimensional case. These criteria only involve the velocity field, which 
indicate that in spite of the nonlinear coupling between the equations for velocity field and 
the phase function, the velocity field indeed plays a dominant role in regularity for solutions 
to system (|1.4p ~ ()1.6p . just like the decoupled incompressible Navier-Stokes equations. 

First, we provide a result on regularity criteria in terms of the velocity u [27] or its gradient 
Vu [2]. 

Theorem 4.1. Suppose n = 3. For {uq,cI)q) € V x Hp, let (n(t),0(t)) be a local smooth 
solution to the problem (|1.4p - (|1.8p on [0,T) for some < T < +oo. Suppose that one of the 
following conditions holds, 

fT 3 2 3 

(i) / \\Vu{t)\\lpdt < +CX), for - + - <2, - <p < +00, 
Jo p s 2 

3 2 

(a) / \\u(t)\\lpdt < +00, for - + -< 1, 3 < p < +00. 
Jo P s 

Then {u{t),(p{t)) can be extended beyond T. 

Proof. We keep in mind that uniform estimates (|2.19p and (j3.7p still hold for t > 0. 
Suppose that (i) is satisfied. For p > |, we estimate ()3.16p as follows 

(w-Vn, An) = — V kUjV jUiV kUidx — / UjV j{V kUi)V kUidx 

Jq Jq 
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Vfc Uj V j UiVkUi dx 



< CIIV-uIIlpIIVu'i^ 



2p 

2p-3 . ,,3 



< CIIV-uIIlp ||Vn|| p ||An||p + ||Vn 



< ^\\Auf + C l^\\Vu\\-LP + \\Vu\\l'f''y\Vu\f. (4.1) 
Combining (|4.1|) with the other estimates in the proof of Lemma 13.31 we obtain that 



d 



6E 



2 / 2p 

I 2p-3 



—A{t) + n\\Au\\^ + kjET] A— <C 1^1 + \\Vu\\ll-' j A{t). (4.2) 

Then by the Gronwall inequahty, we see that A{t) < Ct for t G [0, T], which imphes that 
the X H'^ norm of the strong solution (n, (/>) is bounded on interval [0,T]. This yields that 
[0, T) cannot be the maximal interval of existence, and the solution (n, </>) can be extended 
beyond T. 

Next, we suppose that (ii) is satisfied. For p > 3, we estimate (|3.16p in an alternative 
way such that 

{u-Vu,Au) < C||u||lp||Vu|| 2p IIAiill < ClliillLpllAnll f ||An||t||Vu||^ + ||Vu| 

< ^\\Auf + C (^\\u\\f/ +1^ \\Vuf. (4.3) 

then by the Gronwall inequality, A{t) < Ct for t € [0,T], which again yields our conclusion. 
The proof is complete. □ 

As for the conventional Navier-Stokes equations (see, for instance, [5ll35t[55]). we can 
improve the results in Theorem 14. II and obtain some logarithmical-type regularity criteria for 
our phase-field Navier-Stokes system ()1.4p - p.8p . 



Theorem 4.2. Suppose n = 3. For {uq,(Pq) G (V n H^) x Hp, let {u,(p) be a local smooth 
solution to problem ()1.4p - ()1.8p on [0,T) for some < T < +oo. // one of the following 
conditions holds, 



[ TTTT^TW^W^^^ < +°«' for- + -<2, |<P<6, (4.4) 
Jq 1 + ln(e + ||V-u(t)||Lp) p s 2 

(ii) 

r "j^^l^^^ dt <+oc, for^ + -<l, 3<p<+oo, (4.5) 
Jq 1 + ln(e+ ||u(t)||L°°) p s 

then (n, (p) can be extended beyond T. 

Proof. We recall that uniform estimates ()2.19p and ()3.7p still hold for t >0. 

Case (i). Suppose that (|4.4p is satisfied. Applying A to both sides of equation ()1.4p . 
multiplying the resultant with An and integrating over Q, we get 

-^\\Auf + nWVAuf = (V(u ■S/u),VAu) - ( V ( ^V(^ ) ,VAn) -=1^ + 12, (4.6) 
2 dt \ \0(p J J 
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where 



h < C||VA«||(||Vn||24 + ||n||Loo||An||) 

< C||VAn|| ^||Vn||L4 + ||u|||2 Ikllni 11^^ 

< C||VAu|| (^||VAu||5||Vn||i + ||VAu||t ||Vn||t + ||Vnf 

< ^||VAn||2 + C7(||Vt.||io + ||Vtxf). 

o 

Estimate for the term I2 wUl be postponed. On the other hand, we deduce from equation 
(fOD that 



2di 



d 6E 6E 



k 



6E\ k 



(30^-l)5J(0),A 



-M2{B{cl>)-/3) (dtf 



5E 



JE\ 6k 
M4i.(0)(/(^),Af 



6E 



6E 



= E^- 

The term J( can be estimated in the following way: 

6E 



(4.7) 



4 



< 



< 



-ke'j 
7ke"f 



VA- 



2 / 6E\ 
- k£[VA{u-V(l)),VA ^ 



VA— 



H VA(n • V0) 

7 



7A;e7 



JQ JQ 
+C / {VjViyiUkVk<i)fdx + C / (ViVinfcVjVfc0)2dx 

+C / (VjnfcViViVfc0)2dx + C / (ufcVj ViViVfc</))2dx 
JQ JQ 



'Q 

Jkej 



VA— 



+ J'la + ••■ + •^1/) 



where 



JL + J(, < C||Vn||2^i,3||<A||^2,6 <C(||VAtx||i||V^x||U||Vnf) 



< i||VAuf + C7||Vnf , 



^ + J(e < C||Vu||^oo||0||^3 <C(||VA^z||i||Vn||U||V7 

< ^-\\VAuf + C\\Vu\\\ 

J[, < C||VA?xf||V<^||^oo <Ci||VA^x||2, 

J[f < C\\u\\lo.\\<p\\j,, < C\\u\\i^ ( ^-^ \ 1 
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< C(||n||||VAu|| + ll-ul 



2 M„ „i 

+ C Vn 2 VAu 2 



< -\\V/\uf + C 



5E 



6E 



+ llVtil 



Summing up, we get 



J[< 



Jkej 



VA 



6E 



+ (l + Ci)||VAn|r + C 



5E 



+ 



5E 



+ l|Vnf , 



where Ci is a constant depending on \\uo\\, \\(Po\\h^ and coefficients of the system due to 
estimates ()2.19p and ([37 



The remaining terms J2, Jg can be estimated as for J2, J7 in the proof of Lemma [3.3l 

6E 5E 
by using a similar argument with minor modifications (replacing — — in J2, J7 by A-—). 

d(p ocp 

Consequently, 



i=2 



VA— 



+ C||Vuf + C 



< 



3A;e7 „ . (5^ 2 



VA— +C7||Vn|r + C 



A— 



2 



Set 



Now we turn to estimate l2'- 



h < 



< ^i|VAu||2 + C 



V/^uf + \\V(l)\\lo 



6E 



m 



V- 



4(1 + 



< ^||VAdP + ^ 





' + I|A<^II 


2 




2 


5E 


5^ 


2 ^ <5^ 2 



6E 2 

L3 



VA— 



+ 



(5^ 



Collecting the above estimates together, we deduce that 



d 
di 



iu\\ +rji 



V 



6E 



+ 



VA 



6E 



+ ^||VAti| 





5E 


4 


5E 


2 


( 


66 


+ 


66 


+ 



< c 

< C{A^{t) + A{t)) . 
For the sake of simplicity, we denote 

Q{t) = 



+ ||Vn||^° + llVu 



l+l|Vn(t)||^fi 
l + ln(e+ ||Vu(t)||Lp)' 



For I < p < 6 we infer from ()4.2p that 



d 
di 



2p 



A{t) < C7 (^1 + ||Vn||^r'J ^(t) 

< C,Q{t) [1 + In {e + \\Au{t)\\)] Ait), 



(4.8) 



(4.9) 
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1 



where C^, is a constant depending on ||uo||, II^oIIh^ coefficients of the system. 

Because of (03]), we denote Jq Q{t)dt = M < +00. Fix e G (o 
= to < ti < ... < tjy-i < t]\f = T such that 
ri 



/ Q{t)dt < -, 



V i = 1,2, ...,iV, with iV 



2M 



Then there exist 



+ 1. 



Set to = Iq = Q and t^+i = ij^ = T. It follows from our assumption on the initial data 
that .4(0) < +00. Due to ()2.20p . for each z = 1,2, ...,N, there exists ti € such that 

A{ti) < +00. Moreover, 



Q{t)dt < I Q{t)dt < 2 ^ ^' ^ = 0,iV, 



(4.10) 



Q{t)dt< Q{t)dt+ I Qit)dt<e, for i = 1, 2, iV - 1. (4.11) 



ti+i 



We can prove the required result by an iteration argument from i = to i = N. For i = 0, 
it follows from the Gronwall inequality and (|4.10p that 

A{t) < >l(0)exp Tc* [l + ln(e + sup[o,t]||An(-)|l)j / Q{s)ds 



< ^(O)e'^*' ( e + sup[o,i] \\Au{- 



c. 



We infer from (jM]) and (lil^ that for t £ [0,ti], it holds 



|(l|Anm+.i 



Vt G [0,ti]. 



(t)] <C(e + sup[o,j]||An||(- 



(4.12) 



(4.13) 



V- 



.6E 



Since ||Ati(0)|| and 
grating ()4.13p from to t, we get 

,6E 2 



(0) are bounded due to our assumption on the initial data, inte- 



||Au(t)f +7?i 

< ||An(0)|p+r/i 



V- 



(0) + C7r(e + sup[o,i]||A7x(.)||), Vte [0,ti]. 



Then taking the supremum of both sides for t G [0, ti], we can see that 

,SE. 



sup[o,ii] ( ||An(-)r V— (• 
5E 2 



< ||An(0)||2 + r/i 

< \\AuiO)f + m 



V- 



(0) + Cr e + sup[o,^,]||An(• 



l 



(0) + -sup[o,t,]||At/(-)f + Ct, 



(4.14) 



which indicates that 1 1 An 1 1 and 



are uniformly bounded on [0, ti]. 

Then we can repeat the above argument for i = 1, N such that on each interval [ti, tj+i], 
it holds 



6E,^ ^ 



s^p[u,u+i] iiM-)ir + ^1 V— (•) < iiAn(ti)r + 



V 



.6£ 

6^ 



{t,) + CT, (4.15) 
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where the bound of ||Au(tj)||, (^«) given by the estimates in the previous step on 

As a consequence, we get 

II^IIl°°(0,T;H2) < |1'/'||l°°{0,T;H5) < 

which indicate that [0, T) cannot be the maximal interval of existence, and the solution (n, </>) 
can be extended beyond T. 

Case (ii). We re-estimate the terms Ii and J[j in a different way by using the uniform 
estimates (j2.19p and (|3.7p . The estimate for Ii can be done as follows: 

< C||VAu|j ( ||An|| 2p + lltillLP 

V Lp-2 

< CIlVAiillllVAullillAiill^^tllnllLP 



^1 



_2p_ 



< ^llVAuf + C||n| 
Meanwhile, using Lemma |3. 11 we get 

J[f < C||n||2^||</>||2,, <C||n||2^||</.||^5||<A||H3 



< CdlullilVAull + ||n|| 



V- 



.6E 



+ C 



< ^\\VAuf + C 



V 



.6£ 



+ llVul 



The other terms in /2, J{, Jg are estimated in the same way as in Case (i). Then we deduce 
that 



d_ 

di 



|An||^ + Vi 



V- 



.6E 



+ 



VA 



5E 



+ /i||VAu| 



2p_ 

< C\\u\\l-^ \\Auf + C 



< C ( 1 + ||n||£;=' 



_2p_ 





2 


6E 


2 




+ 


66 








_6E 


\\Auf 


+ ??1 







2 



[1 + ln(e + ||Auf )] + \\Auf + , (4.16) 



1 + ln(e + 

where we have used the Poincare inequality 

+cm 

We infer from ()4.16p and the Gronwall inequality that for all t € [0,T], 

5E "2 



6E 


2 


5E 


2 




2 


5E 


< c 




+ c 




< c 




5(t> 















In 1 + ln e+ ||Au(t)|r + r^i 



< Cln ^1 + ln |^e+ |[Ano||^ + 7/1 
which together with (|4.5p implies that 



it) 



(0) 



2p_ 

LP 



1 + ln(e + ||u(t)||L°° 

'"IIl°°(0,T;H2) < C": II</'IIl°°(0,T;H5) < C*- 



-dt, 



The proof is complete. 



□ 
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5 Stability 



Denote the total energy of the system ()1.4p - ()1.8p by 

m = \\Ht)\?+Em))- 

We recah that £{t) satisfies the basic energy law (|2.16p . which characterizes the dissipative 
nature of the problem (jl.4p - (jl.8p . Inspired by [21] for the liquid crystal system, we can show 
that if the initial datum is regular and the total energy 8{t) cannot drop too much for all 
time, then our problem (jl.4p - (|1.8p admits a unique bounded global strong solution. 

Proposition 5.1. Let n = 3. For any initial data (uq, (po) € V x Hp, there exists a constant 
Eq € (0,1), depending on ||wo||hi;II0o||_h'4 and coefficients of the system such that either 

(1) The problem (|1.4p - (jl.Sp has a unique global strong solution {u, (p) with uniform-in-time 
estimate 

\Ht)\w + mt)\\ii^<c, yt>o, (5.1) 

or (2) there is aT^ £ (0, +oo) such that f (T*) < f (0) — Eq. 

Proof. The proof is based on the higher-order differential inequality ()3.14p (cf. Lemma IS.Sp 

and an argument similar to that in |2HI25j. We only sketch it here for convenience of the 

readers. For any initial data (uo,0o) € V x Hp, let L be a constant such that ||V«o|P + 
5E 2 

——(0) < L. It follows from Lemma l3.ll that |[(/>o||/f4 can be bounded in terms of L and 
ocp 

||(^o||_ff2. Then by Propositions 12.11 and (jS.ip . ||ti(t)|| and ||</'(t)||H3 can be bounded by a 
constant depending L, ||(/)o||j72 and coefficients of the system. Then we can fix the constant 
r] in the definition of A{t) (cf. (|3.13p ) and C* in p.l4p . Consider the ODE problem 

jY{t) = C,[(Y{t)f + Y{t)], y(0)=max{l,?7}L>^(0) (5.2) 

The maximal existence time Tmax of the unique local solution Y(t) is determined by 1^(0) 
and C*. Now we take 

to = ^TmaxiY{0),C^), eo = ^ min {/i, 7} . 

If (ii) is not true, we have £{t) > £{0) — £0 for all t > 0. From the basic energy law (|2.16p . 
we infer that 

/ A{t)dt < / A{t)dt < k£q, with k = max{l, r/} max{^~^, 7^^}. 
J'f Jo 

Hence, there exists a G [^,^0] such that A{t*) < < 1"(0). Take t=K as the initial time 
for (|5.2p with Y{t*) = 1^(0), we infer from the above argument that Y{t) and thus A{t) 
remain bounded at least on [0, ^] C [0, t^ + to] with the same bound as that on [0, to] (since 
the bound for Y{t) is the same). An iteration argument shows that A{t) is bounded for all 
t > 0. The proof is complete. □ 

Corollary 5.1 (Eventual regularity of weak solutions in 3D). When n = 3, let {u^cj)) be a 
global weak solution of the problem (ll.4p - ()1.8p . Then there exists a time Tq G (0, +00) such 
that {u,(l)) becomes a strong solution on [To,+oo). 
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Proof. It follows from ()2.19p . ()2.20p and Lemma l3. II that there exist a time Ti > such that 
||«(ri)||, ||(/)(ri) ||j:^4 are bounded. Taking Ti as the initial time, we can fix L in Proposition 
5.11 and thus Eq. (|2.20p yields that there exists a Tq > Ti such that 



< L, 



(^11 Vn| 



To 



6E 



)dt < eo- 



Taking Tq as the initial time, we can apply the argument for Proposition 15.11 that {u, (J)) will 
be bounded in x H'^ after Tq. □ 

Definition 5.1. We say cp* G Hp is a local minimizer of the elastic energy E{(j)), if there 



exists a 6 > 0, E{(f)*) < E[(j)) for all (p S satisfying 



\h2 < 6. If for all (p G H;, 



E( 



< E[(f)), then (p* is an absolute minimizer. 



Lemma 5.1. Let B be a bounded closed convex subset of Hp. The approximate elastic energy 
E{(f)) admits at least one minimizer (p* G B such that E{(p*) = inf £'(0). 

066 



Proof. Since Eicp) > for all (p B and lim 
minimizing sequence (pn ^ B such that 



E{(p) = +00, E{(p) has a bounded 



EicPr,) 



MEU). 

<t>eB 



(5.3) 



Recalling the definition of E{(p) (|1.3p . we can rewrite E in the following form: 



E{<P) = ^\\A^' + Fi<p) 



with 

F(0) 



V(p ■ V{p^ - (p)dx + 



2e3 



Q 



fdx + -MMicP) - af + -M2{B{cP) - {if. 



Since (pn is bounded in iJ^, there is a subsequence, still denoted by (pn, such that (pn weakly 
converges to a certain function (p* in H^. We infer from the compact Sobolev embedding 
theorem in = 3) that (pn strongly converges to (p* in and H^. It turns out that F{(pn) — > 
F{(p*). Since ||A(/)|p is weakly lower semi-continuous, it follows from (j5.3p that E{(p*) = 
m.iE[(p). Using the elliptic estimate and a bootstrap argument, we see that the minimizer 

(p* is in fact smooth. □ 

Remark 5.1. If (p is a minimizer of E{(p), then it is a critical point of E{(p). It is easy to 
verify that any critical point of E{(p) in Hp is equivalent to a weak solution to the forth-order 
nonlocal elliptic problem 

6E 

— = 0, with (p{x + a) = (p{x) . (5.4) 



In order to prove our stability result, we recall the following Lojasiewicz-Simon type 
inequality whose proof is postponed to the next section. 

Lemma 5.2. Suppose n = 3. Let ip be a critical point of the elastic energy E. There exist 
constants /3 > 0, € (0, ^) depending on ip such that for any <p G Hp{Q) with \\(p — 'ip\\u2 < (3, 
it holds 

>\E{cP)-Em^-''. (5.5) 



5E 
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Now we state the main result of this section. 

Theorem 5.1. Let cp* € Hp(Q) be a local minimizer of E((j)). For any R > 0, consider the 
initial data 

{uo,(Po)^B = {{u,(P)eYxH^{Q): \\u\\ni < R, Uo - (I^*\\h^ < R}- 

For any e > 0, there exists a constant a £ (0, 5) that may depend on (p* , R, e and coefficients 
of the system such that if the initial data {uo,(po) G B satisfies the condition 

IWoW + llf^o - < (5.6) 

then problem ()1.4p - (|1.8p admits a unique global strong solution satisfying 

-<p*\\H2<e, Vt>0. (5.7) 



Proof. If 1 1 uq 1 1 HI ^ R and ||(/)o — (/>*||//4 < i?, then the constant in Proposition 15. II depends 
on (p^:, R and coefficients of the system. It follows from Proposition 12.11 and Proposition 13.11 
that ||n(t)|| and ||'/*(i) ||_f/3 uniformly bounded (by a constant depending on (p*, R and 
coefficients of the system). In what follows we denote by C, Ci genetic constants that only 
depend on R, cp^^ and coefficients of the system. 
By a direction computation, we get 

E{cPo) - E{cP{t)) = [E,{cPo)-E,{cPo)] + ^Ah[{A{cPo)-af-{AicP)-af] 

+^M2[{B{<Po) - f3f - {B{cP) - P^f] 
'■= El + F2 + F3, 

where 

El < C{\\A^o-^n + \\{(po-(p){4 + ^^ + <Po^ + m) 

< C{\\A^o - A<^i| + \\<Po - (PWUl + 0' + M + 

< CUo-Hm, (5.8) 
F2 < C||</.o-0||Li(||0o + 0||Li+2|a|) <C7||0o-0||, (5.9) 
E3 < C||V(/>o + V0||||V0o-V</.||+C|l,/.o + (/>||||(/>g + </.'-2||L.o||(/>o-0|| 

< C\\cPo-(P\\m. (5.10) 

Since the total energy £ is decreasing in time, we infer from the above estimate that 
0<^(0)-^(t) = ^\\uof-^\\u{t)f + E{cPo)-E{cP{t)) 

< l\\uof + E{<Po) - E{(P{t)) 

< l\\uof + Ci\\cP{t)-cPo\\H2. (5.11) 
Let /3 denote the constant in Lemma 15.21 that depends only on ip = (p*. Denote 

t^ = mm|l,eo,e,-,-,— 
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We assume that a < jw. Let T be the smallest finite time for which ||(/>(T') — (j)*\\fj2 > w. 

Then by the proof of Proposition 15.11 the problem admits a bounded strong solution on 

[0,T). If there exists G (0,T) such that = E{(j)*), since cj)* is the local minimizer and 

§E 

Uit*) - 4>*\\m <TU<6,we deduce from ([236]) that ||Vu(i)|| = — (t) 

follows from 



for t > U. It 



k\\ < \\u-V(j)\\+-/ 
< Cf||Vn|| + 



6E 



< CIlV-ullllV^llLa +7 



5E 



SE 



(5.12) 



that for t > t*, \\(l)t\\ = 0. Namely, (p is independent of time for t > t^. As a result, u(t^) = 
and (j){t^) = (p** , where (p** is also a local minimizer (but possibly different from (p*). Due 
to the uniqueness of strong solution, the evolution starting from will be stationary. The 
proof is complete in this case. 

We proceed to work with the case that £{t) > E[(p*) for all t £ [0, T). From the definition 
of T, we see that the conditions in Lemma 15.21 are satisfied with ^p = cp*, on the interval 
[0,T). Consequently, 



/i||Vnp + 7 



- -ig(t) - Eir)f > 



6E 



i||n||2(i-«) + 



6E 



> C \\Vu\\ + 



5E 



, V t G [0,f). (5.13) 



We infer from (|5.12p that 

r Utmdt < C{£{0)-E{(P*)f <C{\\uof^ + \E{(Po)-E{ 
Jo 



< C(||nof' + |l</>o 



which implies that 



\\(p{T) - (P*\\h^ 

< ||(/>(f') - (/>o||/^2 + 1100 - (/'*||//2 

< cum - Mhmf) - <Pof^ + Uo - ^*\\h^ 

1 

cT \ 3 



< c 



k{t)\\dt\ +\\(Pq-(P*\\h2 



< C2 \\uo\\ 3 + \\(po 



+ \m 



Taking 



a < min 



4 ' V4C2 



(5.14) 



(5.15) 



we easily deduce from (|5.14|) that ||0(r) — ^ n^^^ ^ ^ 
with the definition of T. Hence, T = +00 and there holds 



< jzu < vj, which leads to a contradiction 



\u2 < < e, V t > 0. 



(5.16) 
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As a consequence, 

Hit) - Mm < Mt) - + U* - 4>o\\h^ <w + a<^w. (5.17) 

then it follows from ()5.1ip that 

£{t)>£{0)-eo, Vt>0. (5.18) 

By Proposition lS-H we see that ()1.4p - (|1.8p admits a unique global strong solution that satisfies 
()5.7p . The proof is complete. □ 

Corollary 5.2. Assume that the assumptions of Theorem \5.1\ are satisfied. The global strong 
solution {u, (j)) has the following property: 

lim {\Ht)\\^. + \\(t>{t)-<t^oo\\m)=^, (5.19) 

where (j)oo G Hp is a solution to ()5.4p such that E{(f)*) = £'((/>oo). Moreover, there exists a 
positive constant C depending on UQ,(f)Q and coefficients of the system such that 

\Ht)\\m + Uit) -^ooWm <C{l + ty^^, Vt>0. (5.20) 
d' G (0, ^) is the Lojasiewicz exponent in Lemma \5.S\ depending on (f>ao- 



Proof. We infer from the higher-order energy inequality ()3.14p and uniform estimates (|5.ip 
that ■^A{t) is bounded for t > 0. On the other hand, the basic energy law implies that 
A{t) € L^(0,-|-oo), then we have limt_j._|_oo -4(t) = 0. Thus, we obtain the decay property of 
the velocity field u in V and 

hm —it) =0. (5.21) 

t—>-+oo ocp 

Recalling the proof of Theorem 15.11 we have shown that ||(?i't(t)|| € L^(0,-|-oo). As a conse- 
quence, 4){t) will converge in to a function (/)oo € Hp{Q) that satisfies (|5.4p due to ()5.2ip . 
It follows from ()5.16p that for sufficiently large t, we have 



< Uoo-m\\m + \m)-r\\m 



1 



< min{(5, /3}. 



Thus, applying Lemma 15.21 with (p = (j)oo and = (J)*, we have 

= 0. 



The limit function (j)oo is also a local minimizer of the energy E and it will coincide with 
if the latter is isolated. Finally, the proof for convergence rate (j5.20p is based on Lemma [57 
and higher-order differential inequalities as for the liquid crystal systems [251133]. Since the 
proof is lengthy but standard, we omit the details here. □ 
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Remark 5.2. We can also prove the long-time behavior for global weak solutions with ar- 
bitrary initial data in contrast with smallness assumption like (|5.6p . Indeed, Corollaru I5.il 
implies that any global weak solution {u,(p) to the problem (|1.4p - ()1.8p will become a bounded 
strong one after a sufficiently large time. Then we can just make a shift of time and con- 
sider the long-time behavior of bounded strong solutions. Applying Lemma \5.2l we can use 
the Lojasiewicz-Simon technique (cf. e.g., ll[\^ lF/\[2m.[28[ [3 ^ for various applications) to 
show that each weak solution does converge to a single pair (0,(^oo) with 4>oo satisfying the 
stationary problem (|5.4p . Besides, one can obtain the estimate on convergence rate as (|5.2Up . 



6 Appendix: The Loj as iewicz— Simon inequality 

In this section, we prove that a Lojasiewicz-Simon type inequahty holds in a proper neigh- 
borhood of every critical point of the functional E{(f>). 

First, We recall the definition of analyticity on Banach spaces [371 Definition 8.8]: Sup- 
pose X, Y are two Banach spaces. The operator T : D{T) C X — > y is analytic if and only 
if for any xq (z X, there exists a small neighborhood of xq such that 

T{xo + h)-T{xo) = Y,Tnixo){h,...,h), V/iGX, \\h\\x < r « I. 

n>l 

Here T„(xo) is a continuous symmetrical n-linear operator on Y and satisfies 

XI \\'^n{xo)\\L(^xr^,Y)\Mx < +00- 
n>l 

Proposition 6.1. Suppose n = 3, we have 

(1) E{(l)) : H^iQ) is analytic; 

(2) f| : HjXQ) LliQ) is analytic; 

(3) for any G Hp{Q) , E"{iIj) is a Fredholm operator of index zero from Hp{Q) to Lp{Q). 



Proof. (1) It follows from the definition of E that it is the sum of integrations of polynomials 
in terms of A0, Vcp and cp. Since cp ^ Hp{Q), then those functions belongs to Hp{Q) which 
is a Banach algebra for the pointwise multiplication when n = 3. Thus, £'(<^) : Hp(Q) — )• M 
is analytic. 

(2) Recalling ()3.2p . we see that ^ = keA'^(j)-\-H{(j)), where H((j)) is the sum of polynomials 
in terms of Ai;^, V0 and ip that belong to Hp{Q) when (j) £ Hp{Q). Thus, our conclusion 
follows. 

(3) For any ip,wi,W2 € Hp{Q), we calculate that 

E"{'lp){wi,W2) 
= ^{E'illj -\- SWi),W2)\s=0 

as 



d f 

— / {kgiip + swi) + Mi{A{i; + swi) - a) + M2{B{ip + swi) - (3)f{i: + swi)) W2dx 
ds Jq 

f 6k 12k 12k 
/ ke{A'^wi)w2 \V^pfwlW2 {ipVip ■ Vwi)w2 {ipAil))wiW2dx 



s=0 



e e e 
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I — (SV'^ - l){Awi)w2dx + \ [ (ISV-^ - I27p^ - l)wiW2dx + Mi [ widx [ W2dx 

Jq ^ ^ Jq Jq Jq 



+M2 / fWwidx / f{i;)w2dx + M2{B{i;)-P) 
JQ JQ 

ke{A'^wi)w2dx + R{ip){wi,W2)- 



-eAwi + ^(SV'^ - l)wi 



W2dx 



We observe that the leading order term of the hnear operator E"{'ip) is feeA^. This forth-order 
operator (subject to periodic boundary conditions) from H^{Q) to Lp{Q) can be associated 
with the symmetric bilinear form A : Hp{Q) x Hp{Q) M given by 



A{f,g) = ke [ AfAgdx, yf,ge H^{Q) 
JQ 



such that by integration by parts JgheA"^ f gdx = A{f,g) for any f,g & Hp{Q). Obviously, 
A{-, •) is bounded on Hp{Q). By the elliptic estimate (|2.10p . for any A > 0, G Hp{Q), there 
is some r]' such that A{(p,(p) + ?/|I</>|p > ??'||<^||/^2. Thus, it follows from the Lax-Milgram 
theorem that the self-adjoint operator keA + rjl : Hp{Q) — )■ Lp{Q) is an isomorphism. Then 
we see that keA'^ : Hp{Q) — > Lp{Q) is a Fredholm operator of index zero (cf. e.g., [371 Section 
8]). The remaining term R{'ip) consists of tp, VV', Aip and their integrals as well as differential 
operators A, V. Therefore, R{'ip) is a compact operator from Hp(Q) to Lp{Q) for ij; G Hp{Q). 
they rms of A(/), and 4> that belong to i?p(Q) when G H^{Q). 

As a consequence, for any -0 G Hp{Q), E"{il)) is a compact perturbation of a Fredholm 
operator of index zero from Hp{Q) to Lp{Q), then itself is also a Fredholm operator of index 
zero from Hp{Q) to Lp{Q) (cf. e.g, [371 Section 8]). The proof is complete. □ 

Using Proposition 16.11 we can infer from the abstract result [6l Corollary 3.11] that the 
following result holds 

Theorem 6.1. Suppose n = 3. Let tp he a critical point of energy E. Then, there exist 
constants /3i > 0, G (0, |) depending on ip such that for any 4> G Hp{Q) with \\cf)—ip\\uA < j5i, 
there holds 

f)E 

>\E{cP)-Em^-'. (6.1) 

Proof of Lemma 15. 21 Based on Theorem 16. 11 we now relax the smallness condition and 
show that ()5.5p holds if one only requires that cj) falls into a certain ^-neighbourhood of ip. 
For any (p G Hp{Q), using the regularity theory for elliptic problem, we can see that 



<Mi\\A\<p-i^)\\ + \\<P-i;\\), (6.2) 



where M is a constant independent of (p and tp. If ||i?i>— i/;||j:^2 < 1, we take this assumption just 
to ensure that that the fact ||0||iy2 < HV'Hifa + 1 depends only on tp. Similar to (|5.8p - (|5.10p . 
we see that 

\E{cP)-Em^-' <Ci\\cP-i^\\]^-r (6.3) 
By Holder inequality and Sobolev embedding theorem Hp{Q) ^ L°°{Q), we get 



\Bi<p)-Bii^)\ < C||V(0 + V^)||||V(,^-^)|| + C(||</.||ioo + ||V||i^ + 1)110 -VI 
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\\f{4>-fm\ < C\\AcP-A^P\\ + Ci\\^lo. + \ml^+l)U-^P\\<CU-iJ\ 



where C only depend on ||'0||j:/2. Recalling the expression of H((j)) given in ()3.2p . we obtain 

\\H{cP)-Hm 

< CU - V'IIl- ||V<A||i4 + C7||V||l-||V(<A + V')llL4i|V(<A - 
+C||<A||ioo||A(<A - + IIA^IIII,^ + i^h^U - i^h^ 

+C{U\\l^, IIV'IIl-)||0 - V'll + M,U - ijhi + M2\B{4>) - B{^)\\\fm\ 

+M2(|i?(V)|+/3)||/((A)-/(V')ll 

< c2u-nH^. 

Since Ci,C2-the constants above-only depend on ||'i/'||j:^2 , there exists a (sufficiently small) 
constant /3 independent of (j) which satisfies 

\ n Pi ( Pike 



2M' \ AM{Ci + C2) 



such that if Wcj) — "011^2 < /3, then 



WHic^) - nm + \E{^) - Em'-' < (6.4) 

For any (j) € Hp{Q) satisfying ||0 — V'IIh^ < there are only two possibilities: (i) if (p also 
satisfies ||(/) — V'll//* < then ()6.ip holds; (ii) otherwise, if ||(/) — V'IIh* ^ /^i, noticing that ip 
satisfies (15.41). we deduce from (16.211 and (16.41) that 



5E 



= \\keA'{(P - i;) + H{4>) - H{ij)\\ > ke\\A'{4> - - \\H{(P) - Hm 

> ^>m)-Em'-'. 

The proof is complete. 
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